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Abstract 



We consider QCD corrections to Higgs boson production through gluon-gluon fusion in 
hadron collisions. We compute the cross section, performing the all-order resummation 
of multiple soft-gluon emission at next-to-next-to-leading logarithmic level. Known 
fixed-order results (up to next-to-next-to-leading order) are consistently included in 
our calculation. We give phenomenological predictions for Higgs boson production at 
the Tevatron and at the LHC. We estimate the residual theoretical uncertainty from 
perturbative QCD contributions. We also quantify the differences obtained by using 
the presently available sets of parton distributions. 
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1 Introduction 



Higgs boson production is nowadays a topic of central importance in hadron collider physics 0121 • 
The Tevatron will be actively searching for a Higgs boson signal in the near future. The discovery 
of the Higgs boson is also a primary physics goal of the LHC. 

The main Higgs production mechanism at hadron colliders is the gluon fusion process [Hj, 
an essentially strong-interaction process, that has attracted a large amount of theoretical work 
in recent years. Indeed, limits on the Higgs mass will rely upon QCD calculations of the cross 
sections. Conversely, if a Higgs boson is discovered, discrepancies of its measured cross section 
from QCD calculations may signal deviations of the Yukawa couplings from the Standard Model 
(SM) predictions. It is thus important to provide an accurate calculation of the Higgs production 
cross section, together with a reliable estimate of the associated theoretical error. 

The QCD computation of the gluon fusion production cross section was carried out at next-to- 
leading order (NLO) in Refs. [HE] in the heavy-top limit, and in Ref. [B] with the full dependence 
on the top-quark mass. Perturbative corrections at the NLO level were found to be quite large 
(of the order of 100%), thus casting doubts upon the reliability of the perturbative calculation. 

Next-to-next-to-leading order (NNLO) corrections have been computed in the heavy-top limit. 
The virtual contributions were evaluated in Ref. [Zj. The soft contributions were computed in 
Refs. jSJ El- The remaining hard real contributions were included in Ref. ^U], with a semi- 
numerical approach. More recently, fully analytical results for the real contributions have been 
obtained dHOES]. 

Numerically, it is found that NNLO corrections are moderate in size. There is thus good 
hope that the NNLO calculation gives a reasonable estimate of the cross section. Nevertheless, a 
better understanding of the pattern of radiative corrections would prove useful both to assess the 
reliability of the available calculations and to include higher-order corrections. 

In the present work, we include the dominant effect of the uncalculated higher-order terms, 
by exploiting the resummation of soft-gluon emission. The possibility of performing such an 
improvement relies upon the observation that an accurate use of the soft-gluon approximation 
provides the bulk of the NLO term and a reliable estimate of the NNLO effects |S]. Now that the 
NNLO corrections are fully known, and confirm that observation, it makes sense to add the higher- 
order terms that can be obtained in the soft-gluon approximation, in order to give a more precise 
prediction. Since the soft approximation proves reliable for the NLO and NNLO contributions, 
we make the reasonable assumption that it maintains the same reliability for higher-order terms. 

The paper is organized as follows. In Sect.|2]we give our notation for the QCD cross section and 
the fixed-order radiative corrections. In Sect.|H]we present the formalism of soft-gluon resummation 
to all logarithmic orders, and derive the explicit resummation formulae up to the next-to-next-to- 
leading logarithmic (NNLL) level. 

The resummation formalism correctly predicts the structure and the coefficients of the singular 
terms in the exact NLO and NNLO calculations. The quantitative reliability of the soft-gluon 
approximation can be tested by comparing the truncation of the resummation formalism at the 
NLO and NNLO levels with the exact result. This comparison is performed in Sect. |3J Several 
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variants of the soft-gluon approximation are considered there, in order to justify the validity of our 
approach. Furthermore, in Sect. 14.21 the order of magnitude of the soft-gluon effects is studied, 
by simply considering the value of the short-distance partonic cross section in TV-moment space. 

In Sect.0 full numerical predictions for the Higgs production cross section, including soft-gluon 
resummation at NNLL accuracy and the exact NLO and NNLO contributions, are given, both at 
the Tevatron and at the LHC We also perform a study of the remaining theoretical uncertainty. 
Incidentally, we note that the calculation presented here is the first calculation of a QCD cross 
section at such nominal theoretical accuracy, namely NNLO+NNLL accuracy. 

Finally, in Sect.|H]we give our conclusions. 

In Appendix El we derive a simple prescription to evaluate the large- N Mellin moments of the 
soft-gluon contributions at an arbitrary logarithmic accuracy. The method is a generalization of 
the prescription z N ~ l — 1 — > —6(1 — z — Nq/N). In Appendix lB| we present a formal proof of the 
equivalence of two different formulations of the iV-space resummation formulae. This equivalence 
clarifies the distinction between the all-order logarithmic terms in the soft-gluon resummation 
formalism and the large-order perturbative behaviour due to infrared renormalons. In Appendix IU1 
we show how to compute the logarithmic functions that control soft-gluon resummation. In 
Appendix \D\ we check the numerical convergence of the fixed-order soft-gluon expansion, and in 
Appendix |E] we provide the soft-gluon approximation of the N 3 LO contribution to the partonic 
cross section. 

Preliminary results of this work have been presented in Ref. |14j . 



2 Notation and QCD cross section at fixed orders 

We consider the collision of two hadrons hi and h 2 with centre-of-mass energy yfs. The inclusive 
cross section for the production of the SM Higgs boson can be written as 



M 2 H )=S^ j dx! dx 2 f a /hi(xx, fi 2 F ) f b/h2 (x 2 , fi 2 F ) / dz 
nh Jo Jo 



61 z TH 



XiX 2 J 

■ a® z G ab (z; a s (,4), M 2 H /fi 2 R ; M 2 h /^ 2 f ) , (1) 

where Mh is the Higgs boson mass, th = Mfj/s, and fip and hr are factorization and renormaliza- 
tion scales, respectively. The parton densities of the colliding hadrons are denoted by f a /h{ x i Vf) 
and the subscript a labels the type of massless partons (a = g, qf, q~f, with Nf different flavours of 
light quarks). We use parton densities as defined in the MS factorization scheme. 

From Eq. (J^J the cross section a a b for the partonic subprocess ab — > H+X at the centre-of-mass 
energy s = x±x 2 s = M^/z is 



Ml) = I a^M 2 H G ab {z) = ff W z G ab {z) , (2) 



where the term 1/s corresponds to the flux factor and leads to an overall z factor. The Born-level 
cross section a^ ' and the hard coefficient function G ab arise from the phase-space integral of the 
matrix elements squared. 
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The incoming partons a, b couple to the Higgs boson through heavy-quark loops and, therefore, 
and G ab also depend on the masses Mq of the heavy quarks. The Born-level contribution 
is0 



a 



(0) 



G, 



5> 

Q 



4jg 
Mi 



288tt\/2 

where = 1.16639 x 10~ 5 GeV~ 2 is the Fermi constant, and the amplitude Aq is given by 



(3) 



Ao(x) 




— ITT 



X > 1 
X < 1 



(4) 



l - VT-z 

In the following, Mq = M t or M& denotes the on-shell pole mass of the top quark or bottom quark. 

The coefficient function G a b in Eq. (JTJ) is computable in QCD perturbation theory according 
to the expansion 



G ab {z- a s (/4), M H /fi 2 R ; Ml/^p) = al(fj? R 



n=0 



7T 



G^^Mll^MH^) 



G 



ii 



o|Q4) ~m / M 2 , M 2 

7T ^ M> 



(1) 
ab 



JL- 

,2 ' 



«s(/4) 



7T Z 



r (2) 



"I 2~ ' 2~~ 



+ 0(a|) 



where the (scale-independent) LO contribution is 

G^ b \z) = 6 ag 5 bg 5(1 - z) . 
The terms G^j and G^ give the NLO and NNLO contributions, respectively. 



(5) 
(6) 

(7) 



The NLO coefficients G^} are known. Their calculation with the exact dependence on M t (and 
Mb) was performed in Ref. jH], where it was also observed that the NLO Higgs boson cross section 
is well approximated by considering its limit M t 3> Mh jHEJ- Therefore, throughout the paper 
we work in the framework of the large- M t approximation: we consider the case of a single heavy 
quark, the top quark, and Nj = 5 light-quark flavours, and we neglect all the contributions to 
G^ that vanish when Mu/M t — > 0. However, unless otherwise stated, we include in the full 
dependence on M t and M&. At NLO this approximation jHl E3 turns out to be very good when 
Mh < 2M t , and it is still accurate''' to better than 10% when Mh ^ 1 TeV. 

The use of the large- M t expansion considerably simplifies the calculation of the QCD radia- 
tive corrections, since one can exploit the effective-lagrangian approach [TTl ITHl IT3] to embody 
the heavy-quark loop in an effective point-like vertex. The virtual [7j and soft ^01 1201 EU 122] 
contributions (i.e. the contributions that are singular when z — > 1) to the NNLO coefficients 
G ab (z) were independently computed in Refs. jH] and 0. The hard contributions were consid- 
ered in Ref. [TU], by expanding Gf b '(z) in powers of (1 — z) and evaluating the coefficients of the 

^The accuracy of this approximation when Mh & 2M t is not accidental. In fact, as pointed out in Refs. [51 116| 
and discussed below, the main part of the QCD corrections to direct Higgs production is due to parton radiation at 
relatively low transverse momenta. Such radiation is weakly sensitive to the mass of the heavy quark in the loop. 
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expansion explicitly up to order (1 — z) 16 . An independent NNLO calculation, which includes 
all the hard contributions in closed analytic form, was carried out in Ref. [TT]. This result has 
recently been confirmed ^3] by using a different method of calculation. In the high-energy limit, 
the results of Refs. agree with the calculation of Ref. based on fc^-factorization [23]. 

The NLO coefficient functions G^ b in the large- M t limit (i.e. neglecting corrections that vanish 
when M H /M t -> 0) are [HE] 

G$(z; = 5(1 - z) (y + 6C(2) + In j|) + 12 V,{z) 

11 (1 - z) 3 



In z 



G^(z;M 2 H /fi 2 R ;M 2 H /f, 



ZfJLp 



6- 

l-z 2 
2 (1 - ^ 2 



+ 3 Z 



z 



G^{z-Ml/&-Ml/ti) = f 7 ^± 



G^(z;M 2 H /^ 2 R ;M 2 H /^ 2 F ) = 



(8) 

(9) 
(10) 



where £(n) is the Riemann zeta-function (C(2) = 7r 2 /6 = 1.645 . . . , £(3) 



defined 



ln^l - z) 



1.202 . . . ), and we have 
(11) 



The kernels P a b{z) are the LO Altarelli-Parisi splitting functions for real emission, 
P 99 (z) = 6 



-— + -- 2 + z{l-z) 
1 — z z 



W) 3 z 



12) 



and Pgg g (z) is the regular (when z — > 1) part of P gg (z): 
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1 - 2 

(2) 



(13) 



(2), 
aft 1 



The analytic formulae for the NNLO coefficient functions G ab are given in Refs. fUEI] (V, 

f2) ■ -i 

zG ab (z), according to the notation of Ref. [TT]). 

For the purpose of the discussion in the following sections, we note that we can identify three 
kinds of contributions in Eqs. (jKj)- (flT)|) and in the analytic formulae for G ~^\ 



• Soft and virtual corrections, which involve only the gg channel and give rise to the T>i and 
5(1 — z) terms (see Eq. (JHJ)). These are the most singular terms when z — > 1. 

• Purely collinear logarithmic contributions, which are controlled by the regular part of the 
Altarelli-Parisi splitting kernels (see Eqs. (JBJ), ©)• The argument of the collinear logarithm 
corresponds to the maximum value (^ max ~ (1 — z) 2 Mfj/z) of the transverse momentum 
qx of the Higgs boson. These contributions give the next-to-dominant singular terms when 
z -> 1. 

• Hard contributions, which are present in all partonic channels and lead to finite corrections 
in the limit z — > 1 . 
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In this work we are mainly interested in studying the effect of soft-gluon contributions to all 
perturbative orders. Soft-gluon resummation has to be carried out in the Mellin (or iV-moment) 
space We thus introduce our notation in the iV-space. 



We consider the Mellin transform apf(M^) of the hadronic cross section a(s,M^). The N- 
moments with respect to th = M# fs at fixed Mh are thus defined as follows: 

a N (M 2 H ) = [ dr H r^ 1 a(s,M 2 H ) . (14) 
Jo 

In iV-moment space, Eq. (JJJ takes a simple factorized form 

°n-i{M 2 h ) = a® fa/h u N^l) fb/h 2 N{fi 2 F ) G abiN {a s {fi 2 R ),M 2 H /fi 2 R ;M 2 H /fi 2 F ) , (15) 

a,b 

where we have introduced the customary iV-moments of the parton distributions (f a /h,N) an d of 
the hard coefficient function (G^n): 

fa/h,N{V 2 F )= dx X N ~ 1 f a /h(x,fi 2 F ) , (16) 

J 

G ab , N = f dzz N - 1 G ab (z) . (17) 



Jo 

Once these A^-moments are known, the physical cross section in x-space can be obtained by Mellin 
inversion: 

a(s,M 2 H ) = a^ ^ (^l\ N+1 f a/hl ,M) h, h M) 

a,b J CMP-i^ Zm V S J 

x G abyN (a s (fi 2 R ), Ml/ fa M 2 H I 'fa , (18) 

where the constant Cmp that defines the integration contour in the iV-plane is on the right of all 
the possible singularities of the iV-moments. 

Note that the evaluation of G ab (z) in the limit z — > 1 corresponds to the evaluation of its N- 
moments G a b,N in the limit N — > oo. In particular, the soft, virtual and collinear contributions to 
^ab ( z ) ^ ea( i to In iV-enhanced contributions in iV-space according to the following correspondence 
(see Appendix |A"| : 



i 



dz z N - x V k {z) = ^ ^ fc+1 ln fc+1 N + 0(ln k N) , (19) 
dzz N - 1 8{\-z) = \ , (20) 
dzz"- 1 \n k (l-z) = t^\n k N + o(^\n k - l N) . (21) 
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3 Soft-gluon resummation 



3.1 Summation to all logarithmic orders 

In this section we consider the all-order perturbative summation of enhanced threshold (soft and 
virtual) contributions to the partonic cross section for Higgs boson production. The threshold 
region z — > 1 corresponds to the limit iV — > oo in iV-moment space. We are thus interested in 
evaluating the hard coefficient function G a b,N, by keeping all the terms that are not vanishing 
when iV — > oo. To this purpose, we first note that in this limit only the gg partonic channel is 
not suppressed. In other words, we have: 

G ab MM»R),M%/ f il;M%/&) = 0(l/N) (ab^gg) , (22) 

where the notation 0(1/N) means that the right-hand side vanishes at least as a single power of 
1/N (modulo IniV corrections) when N — > oo. When the partonic channel is ab ^ qq, Eq. (|22|) 
simply follows from power counting: the final state X in the partonic subprocess ab — > H + X 
contains at least a fermion, and the corresponding cross section thus vanishes in the soft limit. 
When ab = qq, the threshold limit selects the exclusive subprocess qq — > H that vanishes since 
the gluon-mediated production of a spin particle through qq annihilation is forbidden in the 
massless quark case. As a matter of fact, gluonic interactions conserve helicity, so that the total 
spin projection along the incoming qq direction is ±1, which is incompatible with the production 
of a spin state. 

Observe that the large- N behaviour of the hadronic cross section <ttv(M^) also depends upon 
the large- N behaviour of the parton densities, according to Eq. (|15|). Thus, the 0(1/ N) relative 
suppression of the partonic cross section in the qg channel relative to the gg channel may be 
compensated by the enhancement of the quark with respect to the gluon density. Under the 
typical assumption that the gluon density is softer than the (valence) quark density at large x, it 
is possible to show (see Sect. 2.4 in Ref. |2Zj) that the two parton channels contribute with the 
same power behaviour in N to the total hadronic cross section. In the present work, however, 
we are not considering the large- N limit of the hadronic cross section <jN(Mj I ) for Higgs boson 
production. We are rather using the soft-gluon approximation to find a good approximation of 
the full partonic cross section, to be convoluted with the parton densities and used in kinematical 
regimes* that are far from the hadronic large- iV limit. In this context, Eq. (|22|) implies that the gg 
channel strongly prevails over the other channels in the evaluation of the cross section for Higgs 
boson production. 

We are thus led to consider the gg partonic channel. The formalism to systematically perform 
soft-gluon resummation for hadronic processes, in which a colourless massive particle is produced 
by qq annihilation or gg fusion, was set up in Refs. [213 123 EH| In the case of Higgs boson 
production, we have 

{+oo 2n ^ 

1 + E «s E ^ ln m iV + 0(1/ N) = Gf£> N + 0(1/ N) , (23) 
n=l m=0 ) 

■fMore discussion about this point can be found in Sect.0| 
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where the non-vanishing (singular and constant) contributions in the large-iV limit can be orga- 
nized in the following all-order resummation formula: 

< . 

(24) 



following all-order resummation formula: 

?S(«s(/4), Mh/i& Mil ^l) = a 2 (» 2 R ) C gg (a s (fi 2 R ), M%/ fi%; M 2 H /&) 

■ exp{fe(a s (/4), IniV; Mf //4 M 2 H /fi 2 F )} 



The function C gg (as) contains all the contributions that are constant in the large- N limit. 
They are produced by the hard virtual contributions and non-logarithmic soft corrections, and 
can be computed as a power series expansions in «g: 

+00 , , 2 \ \ n 

C gg (a s (v 2 R ),M 2 H /» 2 R ;M 2 H /ti) = l + Yl ^\M 2 /^M 2 h /^ 2 f ) . (25) 

n=l ^ ' 

All the large logarithmic terms a§ ln m N (with 1 < m < 2n), which are due to soft-gluon 
radiation, are included in the exponential factor expQn- It can be expanded as 

g H (as^ 2 R )MN^,^f) = ^a s ^^ln-iV (26) 

= \nNg%\b a s (fi R ) In N) + g {2 \b a s (fi 2 R ) In TV, Mf//4; Mf//4) 

+ <xb(Pr) IniV, M 2 H /^ 2 R ; M 2 H /ft 2 F ) 

+00 

+ E M/4)p 2 In AT, M 2 H /» 2 R ; M 2 H /^ F ) , (27) 



n=4 



where, for later convenience, we have introduced the first coefficient, b , of the QCD /^-function. 
The functions gffl are defined such that g^\b as IniV) = when a s = 0. 

Note that the exponentiation in Eqs. (|24|) and (J26|) is not trivial |23l l2*o] . The sum over m in 
Eq. (J23j) extends up to m = 2n, while in Eq. (|26|) the maximum value for m is smaller, m < n + 1. 
In particular, this means that all the double logarithmic (DL) terms a^G^' 2 ^ ln 2n N in Eq. (|23|) are 
taken into account by simply exponentiating the lowest-order contribution agG^' 2 ^ In 2 N. Then, 
the exponentiation allows us to define the resummed perturbative expansion in Eq. ([27)1 . The 
function In N g^ resums all the leading logarithmic (LL) contributions «g ln n+1 N, gj/ contains 
the next-to-leading logarithmic (NLL) terms agln n iV, asQii collects the next-to-next-to-leading 
logarithmic (NNLL) terms «g +1 ln n N, and so forth. Note that in the context of soft-gluon resum- 
mation, the parameter ctg In TV is formally considered as being of order unity. Thus, the ratio of 
two successive terms in the expansion (|27|) is formally of 0(ots) (with no IniV enhancement). In 
this respect, the resummed logarithmic expansion in Eq. (J27|) is as systematic as any customary 
fixed-order expansion in powers of as- 

The purpose of the soft-gluon resummation program is to explicitly evaluate the logarithmic 
functions g^ n > of Eq. ()27|) in terms of few coefficients that are perturbatively computable. In the 
case of Higgs boson production, this goal is achieved by showing that the all-order resummation 
formula (|24|) can be recast in the following form jS]: 

g3KW).«;</^) = 4^l)C gg (a s ^ 2 R ),M 2 J^ 2 R ,M 2 H /^ 2 F ) 

■ A^(ag(/4), MHnl) + 0(1/N) . (28) 
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The factor C 99 («g) in Eq. ("2*8*)) is completely analogous to the factor C gg (as) in Eq. ("2*""]). The 

difference between C gg and C gg is simply due to the fact that some constant terms at large N 

J gg 

representation: 



have been moved from C gg to Aft. The Sudakov radiative factor A^ has the following integral 



Af(a s (/4),^;^) = exp<{ / dz 



2 71^2 \ ( rl ~N-1 



z 



WM'h d q 2 



*4 I 2 



^A(a s (q 2 ))+D(a s ((l-zYM 2 H )) 



(29) 



where A(ac$) and D(a§) are perturbative functions 



A(a s ) = ^ (v)" AW = + (?) ^ + (?) A(3) + ° ( ^ } ' (30) 

n=l 

D(a s ) = Y,(-) 0<">=(*) C' 2 ) + 0(al) . (31) 

n=2 

The coefficients A^ n ^ and are perturbatively computable. For example, they can be extracted 
from the calculation of G g9y N at N n LO. 

By inspection of z and q 2 integrations in Eq. ("2*"")) . it is evident that the radiative factor leads to 
the logarithmic structure of Eq. (|27[). plus corrections of 0(1/ N) that vanish when N — > oo. The 
functions depend on the coefficients in Eqs. ("3~""") and (jUJ), and the functional dependence is 
completely specified by Eq. ("""D*)). More precisely (see Eqs. ("3~9~ |) -(""TI ") ). the LL function depends 
on A^ l \ the NLL function g% depends also on A^ 2 \ the NNLL function gff depends also on A^ 
and D^ 2 \ and so forth. In Appendix El we describe in detail a method to obtain the functions 
g^ (for arbitrary values of n) from the integral representation in Eq. (|29jl. 

The structure of Eq. (J29j) is completely analogous to that of the radiative factor of the Drell- 
Yan (DY) process. The derivation of this result up to NLL accuracy (i.e. keeping only the 
coefficients and A^) was discussed in Refs. [23 123 I2H] • The discussion of Refs. 
can be extended to any logarithmic accuracy by taking into account the following two main points. 

First, beyond 0(ag) the Sudakov radiative factor acquires an additional contribution [3T"1 1311 
13*21 13*3*] due to final-state soft partons emitted at large angles with respect to the direction of 
the colliding gluons (or of the colliding qq pair, in the case of the DY process). In Eq. (J29j) this 
contribution^ is included in the term that depends on the function _D(ag((l — z) 2 Mfj)). 

Second, the term proportional to the function A(as{q 2 )) in Eq. (f2*9*|) embodies the effect of soft- 
parton radiation emitted collinearly to the initial-state partons; it therefore depends on both the 
factorization scheme and the factorization scale \xf of the gluon parton distributions in Eq. (|15[). 
The invariance of the hadronic cross section with respect to fip variations thus implies that the 
function A(as) gets higher-order (n > 3) contributions (given by the coefficients A^ in Eq. flUD) 
to compensate the factorization-scale dependence of the parton distributions, as given by the 

§This contribution is denoted by A^ nt ^ in Refs. |""1 130|. 
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Altarelli-Parisi evolution equations 



d fa/h,N{ni 

din. fip 



5^7ab,iv(as(/4)) fb/h,N(^ F ) • (32) 



It is straightforward to show that the function A(as) in Eqs. (j2HJ) and (jHUj) coincides at any per- 
turbative order with the function that controls the large- N behaviour [33] of the gluon anomalous 
dimensions 7 99 .tv in the MS factorization scheme: 

l 99 ,N(a s ) = -A(a s ) ]nN + 0(l) {N -» oo) . (33) 



The Sudakov radiative factor in Eq. (|29|) can also be expressed by using an alternative integral 
representation: 

r y2M « da 2 ~ 11 

2 / -\A(a s (q 2 )) + D(a s (y 2 M 2 H )) 

(34) 

where Nq = e~ lE (je = 0.5772... is the Euler number), and the new perturbative functions 
C gg (as) and D(as) are completely analogous to the functions C gg (as) and D(as) in Eqs. (|2*H|) 
(unlike C 39 , C 99 does not depend on /x 2 ?) and (|3T|) . respectively. Note that the IniV dependence of 
is fully included in the exponential factor on the right-hand side of Eq. (|34j). The representation 
in Eq. (j3j| was first introduced in Ref. [251 up to NLL accuracy. The equivalence between Eq. (J2HJ) 
and (|34|) to any logarithmic accuracy is proved in Appendix [0 where we also derive the formulae 
that explicitly relate the functions A(a§) and -D(a s ) to the functions C gg (as) and D(as). Using 
Eq. ()34|) . it is straightforward to carry out the z and q 2 integrations and to obtain the logarithmic 
functions gffl in Eq. (|2Tj) (see Sect. 13.21 and Appendix • 

These results on soft-gluon resummation at any logarithmic accuracy deserve some comments 
on the large-order perturbative behaviour. As is well known (see [33] and references therein), any 
perturbative QCD expansion, such as Eq. (J3J), has to be regarded as an asymptotic rather than 
a convergent series. In fact, the perturbative coefficients (e.g. G^) are expected to diverge as 
n\ when the perturbative order n becomes very large. The ambiguities related to the definition 
of the asymptotic series are then interpreted as perturbative evidence of non-perturbative power 
corrections. These features apply to the fixed-order perturbative expansion in Eq. (jSJ) as well as to 
the resummed logarithmic expansion in Eq. (|2*7j) . In other words, the functions gff are expected 
to diverge as n\ at large n. 

Infrared renormalons [33] are a known source of factorially divergent terms in perturbative 
QCD. They arise from the behaviour of the running coupling «s(<7 2 ) in the infrared region and, 
more precisely, from integrating the coupling down to momentum scales q below the Landau pole, 
as set by the QCD scale Aqcd- The representation in Eq. (|2*9*|) may lead to infrared renormalons, 
since it involves z and q 2 integrals that extend in the infrared region independently of the value 
of N. In Ref. jHH], h was observed that the ensuing factorially divergent behaviour corresponds 
to a power correction Aqcd/Mh, which is linear in 1/M#. This observation was based on the 
truncation of the perturbative function A(as) at a fixed perturbative order. However, as pointed 



V /4 »f J 1 Jno/n y 
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out by Beneke and Braun [37|, this simplifying assumption is not sufficient to draw conclusions 
on power corrections. Factorial divergences can arise both from the z and q 2 integrals and from 
the large-order behaviour of the functions A(a^) and D(as): both effects have to be taken into 
account [HZ|. This feature is evident by comparing the integral representations in Eqs. ()29)1 and 
(134)1 . The two all-order representations are fully equivalent, but the integrals in Eq. (|34J) are 
perfectly convergent as long as N < N Mh / Aqcd ~ Mh/ Aqcd- In Eq. (JBH) factorial divergences 
can appear only through the large-order behaviour of the functions A(as) and -D(as) (see also 
Appendix lBj) . The detailed study of Ref. [HZ] shows that the functions D(as) and D(as) (due 
to large-angle soft-gluon radiation) are indeed factorially divergent. In particular, D(as) has a 
factorial divergence that corresponds to a linear power correction and cancels the linear power 
correction found in Ref. Renormalon calculations, based on the explicit evaluation of the 

dominant terms at large Nf, lead to power corrections of the type Aq CD /M h jHZj (or, more 
generally, integer powers of N 2 Aq CD /M h [SHI)- 

Note that the z integral in Eq. (jS3| is not regular when iV > N Mh/ Aqcd ~ Mh/ Aqcd, but 
this does not lead to factorial divergences [39J. We postpone further comments on this point to 
Sect. O 

The all-order resummation formulae presented in this section are valid for Higgs boson pro- 
duction, independently of the use of the large- M t approximation. In particular, the exponential 
factor exp Gh in Eq. (J24*|) and the Sudakov radiative factor Aft in Eqs. (J2HJ) and (JH3J) do not 
depend on M t . The full dependence on M t of the resummation formula (|24)l is embodied in the 
N- independent function C gg (as), namely, in its perturbative coefficients C$ (see Eq. (J25J))- As 
shown below in Eqs. (|44 )) -()47) h in the large- M t limit, the coefficient Cgg becomes independent of 
M t , while Cg 2 , depends logarithmically on M H /M t . 



3.2 Soft-gluon resummation at NNLL accuracy 

In the following we are interested in a quantitative study of soft-gluon resummation effects up to 
NNLL accuracy. We thus need the Higgs bosons coefficients A {1 \ A {2 \ A {3) and £> (2) in Eqs. (JSOI) 
and (|31|). Note that these coefficients are related^ to the analogous coefficients of the DY process 
[23 EDI E] by a simple overall factor C a , which is proportional to the colour charges of the colliding 
partons in the different processes. Therefore we explicitly introduce in the following expressions the 
factor C a , where C a = Ca = N c = 3 in Higgs boson production and C a = Cp = (N 2 — 1)/2N C = 4/3 
in DY production. 

The LL and NLL coefficients and A (2) are well known [321 EES] : 

A« = C a , A® = \c a K, (35) 

with 

™ (!-£»• <»> 

^This relation follows from the general structure |2(JI \22\ of the soft-gluon factorization formulae at 0(a|). 
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The NNLL coefficient was evaluated in Refs. 0E|: 

101 



£ (2) = C a 



C A 



27 



yC(2) + ^C(3) 



14 
27 



C(2) 



(37) 



The higher-order coefficients with n > 3 are not fully known. However, the contribution to 
of the term proportional to NJ~ can be extracted from calculations [SZ1 El] in the large- Nf 
limit. Moreover, by exploiting the relation (j33|) between the anomalous dimensions and A(as), 
the available approximation [33] of the NNLO anomalous dimensions can be used to obtain a 



corresponding numerical estimate 
A® =C a ((13.81 ±0.13) -\ 



of the NNLL coefficient A®: 



55 
48 



-C(3) +C A 



209 5vr^ 



216 54 



+ 



7C(3) \ 



— iV 2 
108 J 



(3f 



13.81 ±0.13) - 2.1467... Nf - 



108' 



-AT 



where we have used the recent analytical computation jJOJ 07] of the term proportional to Nf, 
which agrees with the approximate numerical calculation in Ref. 

The LL, NLL and NNLL functions <7#\<7^ anc ^ 9^h m ^7|) have the following explicit 
expressions (see Appendix : 



*S } (A) 

g%\\,M 2 H liihMlliil) 



+ 



7ib X 
A& 



[2A + (1 - 2A) ln(l - 2A)] , 
2AW lE 



(39) 



[2A + ln(l - 2A)] 



Tlbn 



ln(l - 2A) 



2A + ln(l - 2A) + - ln 2 (l - 2A) 



+ ^-r- [2A + ln(l - 2A)] In 



7r6 



/4 



21^ x , M 
A In 



H 



/4 



(40) 



44(1) 



A 



7T 



2l (1) 7 g&i 
1-2A tt& 2 (1-2A) 



[2A + ln(l - 2A)] 



A^b\ 



+ 



+ 



7T^(l - 2A) 

AW 6a 



2A 2 + 2A ln(l - 2A) + - ln 2 (l - 2A) 



2A + ln(l - 2A) - 

A lE A^ A A^h 



2A 2 



1 - 2A 
1 



2A® A s 



D® A 



ir 3 b 2 1 - 2A n 2 b 1 - 2A 



vr 2 6 
2A^ 



2A 



n 2 b 3 1-2A 



[2A + ln(l - 2A) + 2A 2 ] 



A In 



M 



H 



7T 2 b 

1 /A«6i 



A In 



M 



7T 



/4 



+ 



1 



>A V vr6 2 



[2A + ln(l -2A)] 



2A« 

7T 

4A« 



-A In 



M 



/4 



In 



M 



7T 



7T 2 6 / /4 



+ 



2A« A 2 



7T 



2A 



In 



/4 



(41) 
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where 



■ 2 Mn7V 



A = b a s (/J, R , 

and bo, b±, b<i are the first three coefficients of the QCD /3-function 
1 „ , . 1 



12vr 
1 

W) 



[UC A -2Nf) 
'2857^3 _ 



bi 



24tt 2 



54 



54 



(17C1 - 5C A N f - 3CVA7) 

79 
54 



ClN 



J 



C A N} + y C F iV 2 



(42) 



(43) 



The functions g^' and are well known (see e.g. Ref. [30]). The NNLL function was first 
evaluated in Ref. jlj . Our result in Eq. (|41|) is obtained by using a different method, and confirms 
the result of Ref. |3T| . 



To fully exploit the content of the resummation formula (|24|) up to NLL (and NNLL) accuracy 



we need the constant coefficient C^jj (and Cgg) in Eq. (J23J). The coefficients Cgg and Cgg j read 



,(2) 



(1) 



,(2) 



7T 



67 E ln 



133 7T 2 



+ 



101 

IT 

5N f 7T 2 



14 



M| 
63 



C(3) 



7l 



133 



5N f 



21 vr 2 



7S 11 



2N f 



(44) 



18 T | 



12 



+ 



29 7T 4 



20 



22C(3) - 



4iV/C(3) 



+ ln 2 



M 



H 



/4 



165 



7E + 18 7 | + -N nE + 3tt ; 



3 M 2 M 2 1 

+ £ 7js (33 - 2N f ) In In - -(33 - 2JV,) (67 

l flp flji 4 



+ ln H 



/4 



-36 7 I + (33-2iV / ) 7 | + 7i; 



133 5iV 



7T 2 In 

21 7T 2 



11 7T 2 



N f n 2 



72C(3) 



(45) 



where 



5GW = — 
55 2 



6C(2) 



33 - 2iV 



6 



2 ' 



5G (2) 
gg 



+ h? c(2) _ i c(2)2 _ «» (( 3) 

144 2 sv ; 20 SV ' 4 sv ; 



+ 



19 2 *r 

— + -N f 
8 3 ; 



In 



M 



H 



Mf 



+ N f 



(33 



2iV ' )2 .„4 



48 

169 171 



C(3) 



18C(2)ln 
jj, R n F 

19 



iH£_5 2 £ 3 

144 3 SV 7 6 SV ; 
2 M 2 H 



6 



A7- + (33 - 2JV» C(2) In 



M% 
/4 



465 13 



-JV> - - (33 - 2N S ) C(2) In 



M 
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The terms 5G^j and SG^J are the coefficients of the contribution proportional to 5(1 — z) in the 



.(2) 



(46) 



(47) 



coefficient functions G^J(z) and Gg Z g '(z), respectively. The NLO term SGgg can be read from 
Eq. ©. The NNLO term 8G^ was computed in Refs. 010]. 



.(2) 



(1) 
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As pointed out in Ref. [T^], the dominant part of the corrections of 0(1/ N) in Eq. (J23)) is 
due to collinear-parton radiation, since it produces terms that are enhanced by powers of In N 
(or ln(l — z), as shown in Eqs. (JSJ) and ()21|l). These corrections can be included in the soft-gluon 
resummation formula (|24|1 . In particular, by implementing the simple modification 

OS - OS + MO ^ , (48) 

to the coefficient Cgg , Eq. (|2"4*)l correctly resums O El all the leading collinear contributions to 
Ggg t N, i.e. all the terms of the type (ag ln 2n_1 N)/N that appear in the large- A" behaviour of 
G g n g N . In the following sections we use the prescription in Eq. (f4*K|) to quantitatively estimate the 
dominant corrections to soft-gluon resummation. 



4 Soft-virtual approximation 

4.1 Soft-virtual approximation at NNLO 

In Sect. E] we have discussed a method to resum soft-gluon effects to all perturbative orders. 
The resummed formula is formally justified in the threshold limit z — > 1, where the expansion 
parameter In A" is really large, so that terms that are suppressed by powers of 1/N are negligible. 
We wish, however, to use the resummed formulae also away from the threshold region. Our 
justification for doing so is that we expect that the large- ./V approximation is a good quantitative 
approximation to the exact result. 

The reasons for this expectation are discussed in detail in Refs. [H1E3- We briefly summarize 
the main point of that discussion. In the evaluation of the hadronic cross section in Eq. (|TJ, 
the partonic cross section a a b(s, M H ) has to be weighted (convoluted) with the parton densities. 
Owing to the strong suppression of the parton densities f a /h(x, Hf) at large x, the partonic centre- 
of-mass energy \fl> is typically substantially smaller than y^s ((s) = (xiX2-s) = (r#)s), and the 
dominant values of the variable z = M H /s in the hard coefficient function G a b(z) can be close to 
unity also when yfs is not very close to Mh jlH)- 

At fixed Mh, the reliability of the large- ,2 (large- iV) approximation of the hard coefficient 
function depends on the value of \fs and on the actual value of the parton densities. Therefore, 
the qualitative expectation of the relevance of the large- A" approximation has to be quantitatively 
tested at the level at which the exact result is known, that is up to the NNLO level. Since in the 
rest of the paper we are interested in studying higher-order soft-gluon effects within the iV-space 
resummation formalism of Sect. El in this subsection we study the soft approximations to the 
fixed-order perturbative coefficients in iV-space up to NNLO. At the end of the subsection, we 
also discuss the soft approximations in x-space, which were already considered in Refs. |51 ITBj. 

We begin by defining the A^-space soft-virtual (SV-A r ) approximation of the hard coefficient 
function. We write 



i + G « *>-» + ( ^M>) G msv-N + O(os , 



7T aM ' \ 7T 



(49) 
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The coefficients G^ 2 } T N can be obtained either by Mellin transformation of Gg~g 2 \z), neglecting 



gg,N 

all terms formally suppressed by powers of N, or by expanding Eq. 
We get 



to the fourth order in a$- 



Ml 



G^-'^Ml/^Ml/ti) = Q\n 2 N +\2 lE \nN - 6 lniVln -£ + C<£ , 

fJ> F 



(50) 



GfJ- N {M 2 H /nl, Ml/ti) = 18 In 4 N + In 3 N 
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(33 - 2^)71 + 12 7 f C{? + 7f (g7 ~ y N f - 3 tt 2 ) + H _ 1^ - y C(3) 



.I(33-2JV»ln a ^ 
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i M H 

Mf 



67 5 Ar 3 2 
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+ Ci 2 J . 
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The iV-space s oft-virtual- collinear (SVC-iV) approximation is defined as 



(51) 



G„„ „ [M h /[1r, M H /fi F ) = G N {M H /^ R ; M H /n F ) + 6- 



T 99,N 



GfJ C - N {Ml/&- Ml/ti) = Gfl* Y - N (M 2 H /» 2 R] Ml/ti) + 36 



N 1 
ln 3 iV 
N 



ln 2 iV ln 2 ^, M 2 
72 7i; — 36^^- In ■ 



N 



N 



/4 



(52) 



(53) 



Note that we keep terms proportional to In 2 N/N in the two-loop coefficient G 



(2) SVC-iV 
gg,N 



. These 



subleading collinear terms" appear in the expansion of the (modified) resummed formulae when 
the leading collinear terms are taken into account according to Eq. 



We now want to compare the SV-iV and SVC-iV approximations to the exact results, both at 
NLO and NNLO. To do this, at each order we define the quantity 



A A (fi F , n R ) 



a(n F = n R = M H ) 



(54) 



where the subscript A stands for SV, SVC, or nothing (in the case of no approximation) at the 
given order. The approximated cross sections o"sv an d 0"svc at NLO are defined by Eq. (fTSj) with 
the replacements 



G a b,N 

Ggg,N 



for ah ^ gg 



(4(n 



a s q(i)A-N 



71 



gg,N 



.(2) 



"We recall that these terms are not the complete contribution proportional to In N/N in G gg N . 
anticipate that the effect of including these terms in Eq. l|53|l is numerically negligible. 



(55) 



We also 
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(l) A-N 
gg,N 



where A stands for SV or SVC, and G 

NNLO level, the approximated cross sections a a are defined by the replacements 



is given in Eq. (|5Uj) or (J52*j) . respectively. At the 



G a b,N 
G 99,N 



a s (/x 



R) q(1) 



7T 



ab,N 



for ab 7^ gg 



(56) 



a s(/4) 



i , ^s04) r (i) f «s04) \ 2 ^(2) 
1 + — F~ G «^ + r~^>/ 



A-iV 
N 



where G^^ N is given in Eqs. (JoTj) and (J53*|) . Note that at NNLO, o"sv and o"svc include the 
complete (i.e. without any large-iV approximation and including the qg and qq channels) hard 
coefficient function up to NLO. 

The quantity A is plotted in Figs. ^ and El for the LHC and for the Tevatron, respectively. 
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Figure 1: The SV (dotted lines) and SVC (dashed lines) approximations in N-space versus the 
exact results (solid lines) at the LHC (^/s = 14 TeV). 



The central curves are obtained by fixing ixf = IXr = Mjj. The bands are obtained by varying 
lif and llr simultaneously and independently in the range 0.5M# < ixf,/xr < 2M# with the 
constraint 0.5 < Hf/hr < 2. 

Here and in the following, we use the MRST2002 [5U] set of parton distributions, which includes 
(approximated [IHJ EHl EI|) NNLO parton distributions. The parton densities and QCD coupling 
are evaluated at each corresponding order, by using 1-loop as at LO, 2-loop «s at NLO and 3-loop 
as at NNLO. The corresponding values of as(Mz) are 0.130, 0.1197, 0.1154, at 1-loop, 2-loop and 
3-loop order, respectively. 

As can be observed from Figs. Q and El the SV and SVC approximations in iV-space agree 
very well with the exact NLO and NNLO calculations. Moreover, the differences between the 
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Figure 2: The SV (dotted lines) and SVC (dashed lines) approximations in N-space versus the 
exact results (solid lines) at the Tevatron (yfs = 1.96 TeV). 

approximated and exact results are substantially smaller than the effects produced by the scale 
variations in the exact results at each fixed order. A sizeable part of these small differences can be 
understood from the fact that in the approximation only the gg channel contribution is included, 
thus neglecting the negative contribution from the qg channel, which incidentally is also well 
approximated by its leading collinear behaviour at large N. The very small numerical difference 
between the SV and SVC approximations is consistent with the fact that they formally differ by 
1/N suppressed terms, thus again confirming the consistency of the approximation based on the 
large- N expansion. 

The soft approximation, as well as the SV and SVC approximations, can also be defined in x 
space. In this formulation, the large logarithm to be considered is ln(l — x), instead of In N. It is 
easy to go from one formulation to the other by Mellin transform. For example (see e.g. Eqs. (fT9"j) - 
(|21|) and Appendix El), one goes from the x-space formulae to the iV-space ones by performing a 
Mellin transformation and discarding all 1/N suppressed terms, and all terms that are subleading 
(in the large-iV limit) with respect to the logarithmic accuracy of the initial x-space formulae. 
Thus, x-space and iV-space formulae generally differ by terms that are formally subleading. 

It has been shown in Ref. [39J that large subleading terms arise in the x-space formulation of 
the resummation program. These subleading terms grow factorially with the order of the pertur- 
bative expansion. Their factorial growth is not related to renormalons or Landau singularities. 
These terms are rather an artefact of the x-space approximation, which mistreats the kinemati- 
cal constraint of energy conservation, and they should not be present in the exact theory. As a 
consequence of the factorial growth of these terms, all-order resummation cannot systematically 
be defined (implemented) in x-space, since the series of LL, NLL, NNLL, . . . terms are separately 
divergent. It is therefore interesting to compare the x-space and iV-space approximations for the 
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first few exactly known orders. 
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Figure 3: The SV-x (dotted lines) and SVC-x (dashed lines) approximations versus the exact 
results (solid lines) at the LHC (yfs = 14 TeV). 

At NLO, the x-space soft- virtual contribution (called SV-x approximation, here) to the gluon 
coefficient function G g l g in Eq. (jHJ is 

G$ SY - X (z; Mil fa MH fa = 5(1 ~z)(j + 6C(2) + In 

M 2 

+ 6V (z) \n—f + 12Di(z) . (57) 

The same approximation can be defined at NNLO, and the corresponding contribution to the 
gluon coefficient function G g 2 g(z) has the form 

Gg SV - X (^; MI/aiI, MI/aiI) = 5(1 - z) 5C7( 2 ) + P G< 2) + V 1 Gf + P 2 Gf + X? 3 Gf , (58) 

where the coefficients 8G^ 2 \Gq,Gi and Gg 2 ^ were computed in Refs. [1112] (these coefficients 
can be found in Eq. (2.13) of Ref. |gj). 

After including the leading collinear logarithmic contributions, the x-space soft-virtual-collinear 
(SVC-x) approximation is defined [Hj by 

GCOSVG*^. M 2 7/i 2. ^2) = G (1)SV-X (2 . M 2/^. M 2 //4) _ 12 ln(1 _ z) {5Q) 

G f)WC- X{z . M 2 H/fl 2 R . M 2 h//i 2 f) = G WSV-x {z . M^/ fa, Ml/ fa) - 72 ln 3 (l - z) . (60) 
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Figure 4: The SV-x (dotted lines) and SVC-x (dashed lines) approximations versus the exact 
results (solid lines) at the Tevatron (^/s = 1.96 TeV). 

Figures El and |3] are obtained in the same way as Figs.^and|2l except for the use of the SV-x and 
SVC-x approximations instead of the analogous iV-space approximations. These figures show that 
the bulk of the fixed-order radiative corrections is given by the SV-x and SVC-x approximations, 
as observed in Refs. [H1E3- However, comparing Figs. 01 and H] with Figs. ^ and HI we also see 
that the x-space approximations are worse than the iV-space ones. Furthermore, the difference 
between the SV-x and SVC-x approximations (which is formally suppressed by a power of (1 —z)) 
is considerably larger than the difference between the SV-iV and SVC-iV approximations (which 
is formally suppressed by a power of 1/N). 

The lower quality of the x-space approximations at NLO and NNLO can have two different 
origins. First, the typical value (1 — 2) of the distance from the partonic threshold, which is 
the parameter that formally controls the large-x expansion, can be quantitatively larger than the 
typical value (l/N) of the analogous expansion parameter in the iV-space approach. Second, the 
numerical coefficients in the x-space expansion formulae can be larger than those in the iV-space 
expansion formulae, as is the case at very high perturbative orders, because of the presence of 
factorially-growing subleading terms in the x-space approach. Of course, given the few orders 
available in the perturbative expansion, it is very difficult to explicitly check for the presence or 
absence of these factorial terms. We conclude, however, that these results provide no justification 
for estimating terms of still higher order through the soft-gluon approximation in the x-space 
approach. 
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4.2 Numerical relevance of the resummation 

Before moving to the phenomenological results for soft-gluon resummation in Higgs production 
at the LHC and at the Tevatron, it is interesting to study the main features of the resummed 
coefficient function. This is better done in iV-space, since both the resummed and fixed-order 
coefficient functions in x-space are distributions, and their numerical impact is therefore more 
difficult to assess. 

To directly observe the effect of the exponentiation, we fix the coupling constant at as = 0.1, 
the Higgs mass at Mh = 150 GeV, and the scales at = [j,r — Mh- Figure El shows the 
TV-moments, G g9t N, of the resummed (left-hand side) and fixed-order (right-hand side) gluon 
coefficient function. The fixed-order coefficient function (see Eq. JBJ)) is evaluated at LO, NLO 
and NNLO. The resummed coefficient function (see Eq. (J24)) ) is evaluated at LL (i.e. including 
the function g^), NLL (i.e. including also the function and the coefficient Cgg) and NNLL 
(i.e. including also the function g^ and the coefficient Cgg) order. At large values of N, there is 
a noticeable improvement in the convergence of the perturbative expansion once the resummation 
is performed. Indeed, the ratio between NNLL and NLL results is considerably smaller than the 
one between NLL and LL results. On the contrary, fixed-order results show an increasing ratio 
between two successive orders, which is due to the appearance of new and large ln^ N terms in the 
fixed-order expansion. Note, however, that we do not anticipate very large resummation effects 
on Higgs boson production at the LHC and the Tevatron. In fact, we know [B] fS] that the ratio 
between the NNLO and LO cross sections does not exceed a factor of about 3 at these hadron 
colliders. From inspection of the right-hand side of Fig. we thus expect that the Higgs boson 
cross section at LHC and Tevatron energies is mainly sensitive to the gluon coefficient function at 
moderate values of N. 
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N N 

Figure 5: N -dependence of resummed (left-hand side) and fixed-order (right-hand side) gluon 
coefficient functions for Higgs production (Mh = 150 GeV) with fixed coupling constant as = 0.1. 
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5 Phenomenological results 



5.1 Resummed cross section 



We use soft-gluon resummation in iV-space at the parton level (i.e. at the level of the partonic 
coefficient function G a b) to introduce an improved (resummed) hadronic cross section a^ cs \s,M 2 H ), 
which is obtained by inverse Mellin transformation (see Eq. (|18|) ) as follows: 

r C MP +ioo i N /M 2 \~ N+1 
JCmp-ioo V s / 



X 



Cmp -«oo 



(f.o.) 

+ a^\s,M 2 H ) , (61) 

where a^ Lo '\s, Mfj) is the Higgs boson hadronic cross section at a given fixed order (f.o. = LO, 
NLO, NNLO), is given in Eq. (j2U), and ( G^^ N ) represents its perturbative truncation 

V / (f.o.) 

at the same fixed order in as(/J/ji)- Thus, because of the subtraction in the square bracket on the 
right-hand side, Eq. (joTj) exactly reproduces the fixed-order results and resums soft-gluon effects 
beyond those fixed orders up to a certain logarithmic accuracy. 

In the following subsections, we present numerical results for the resummed cross section 
( j( res )(s ) Mfj) at LL, NLL and NNLL accuracy. The resummed coefficient function G g r ^ N in Eq. (JoTj) 

is evaluated from the expressions in Eqs. (f2"^) - (|2*Tj) : at LL accuracy we include the function g^; 
at NLL accuracy we include also the function and the coefficient Cgg\ at NNLL accuracy we 
include also g^ and Cf g \ Although they are briefly denoted as N fc LL (k = 0, 1, 2), the resummed 
results are always matched to the corresponding fixed order (f.o. = N fc LO) according to Eq. (|61[). 
i.e. LL is matched to LO, NLL to NLO and NNLL to NNLO. 

Unless otherwise stated, cross sections are computed using sets of parton distributions, with 
densities and QCD coupling evaluated at each corresponding order, by using 1-loop as at LO 
(LL), 2-loop a s at NLO (NLL), and 3-loop a s at NNLO (NNLL). 

We recall that the hard coefficient function G a b is evaluated in the large-M t approximation, 
whereas the exact dependence on the masses M t and M& of the top and bottom quark is included 
in the Born-level cross section a^°\ We use M t = 176 GeV and Mb = 4.75 GeV. 

The inverse Mellin transformation in Eq. ()61|) involves an integral in the complex N plane. 
When the iV-moments G gg ^N are evaluated at a fixed perturbative order, they are analytic func- 
tions in a right half-plane of the complex variable N. In this case, the constant Cmp that defines 
the integration contour has to be chosen in this half-plane, that is, on the right of all the possible 
singularities of the iV-moments. 

When Ggg^N is evaluated in resummed perturbation theory, the resummed functions g^\\) in 
Eqs. (JH§j) - (}4*T]) are singular at the point A = 1/2, which corresponds to N = Nl = exp(l/2fc as(/ifj)) 
(i.e. Nl ~ Mh/^-qcd)- These singularities, which are related to the divergent behaviour of the 
perturbative running coupling as near the Landau pole, signal the onset of non-perturbative phe- 
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nomena at very large values of N or, equivalently, in the region very close to threshold. We 
deal with these singularities by using the Minimal Prescription introduced in Ref. In the 

evaluation of the inverse Mellin transformation in Eq. (|61j). the constant Cmp is chosen in such a 
way that all singularities in the integrand are to the left of the integration contour, except for the 
Landau singularity at N = Nl, that should lie to the far right. The results obtained by using this 
prescription converge asymptotically to the perturbative series* and do not introduce (unjustified) 
power corrections^ of non-perturbative origin. These corrections are certainly present in physical 
cross sections, but their effect is not expected to be sizeable, as long as M H is sufficiently per- 
turbative and th = M H fs is sufficiently far from the hadronic threshold, as is the case in Higgs 
boson production at the Tevatron and the LHC. 

The resummed cross section in Eq. can equivalently be rewritten as 

a (res) (s, M H ) = a (sv) 0, M H ) + a ^ ch -) ( s , M%) , (62) 

where a^ SY ^ denotes the contribution obtained by Mellin inversion of G^ es ^, while the matching 
contribution <j( match -) denotes the fixed-order cross section minus the corresponding fixed-order 
truncation of the soft-gluon resummed terms. As can easily be argued from the numerical results 
in Sect. 14. l\ cr^ sv - ) gives the bulk of the QCD radiative corrections to the Higgs boson cross section 
at the Tevatron and the LHC. The order of magnitude of the relative contribution from <j( match -) 
can be estimated from the size of the ratio A (see the definition in Eq. (154)1 ) in Figs. ^and|2J it 
is of 0(10%) and of (9(1%) at NLO and NNLO, respectively. Thus, the fixed-order cross section 
a (match.) quantitatively behaves as naively expected from a power series expansion whose expansion 
parameter is as ~ 0.1. We expect that the currently unknown (beyond NNLO) corrections to 
a (match.) ng^yg no practical quantitative impact on the QCD predictions for Higgs boson production 
at the Tevatron and the LHC. 

We note that the predictions we are going to present regard the production of an on-shell 
Higgs boson. Therefore they are directly applicable at low values of Mh, where the small- width 
approximation is valid. At high values of Mh, corrections due to finite- width effects have to be 
implemented. 

5.2 LHC 

In this subsection we study the phenomenological impact of soft-gluon resummation on the pro- 
duction of the SM Higgs boson at the LHC. The cross sections are computed using the MRST2002 
set of parton distributions |5U|, with densities and QCD coupling evaluated at each corresponding 
order, as stated in Sect. 15.11 and done in Sect. 14. H 

We begin the presentation of our results by showing in Fig. |H]the scale dependence of the cross 
section for the production of a Higgs boson with Mh = 115 GeV. The scale dependence is analysed 
by varying the factorization and renormalization scales around the default value Mh- The plot on 
the left corresponds to the simultaneous variation of both scales, \ip = hr = \Mh, whereas the 
plot in the centre (right) corresponds to the variation of the factorization (renormalization) scale 
1*>f = Xf M h (ji R = xr Mh) by fixing the other scale at the default value M H - 

*An explicit check of the numerical convergence is presented in Appendix iDl 

^The only remaining asymptotic ambiguity is more suppressed than any power law |39| . 
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Figure 6: Scale dependence of the Higgs production cross section at the LHC for Mh = 115 GeV 
at a) (upper) LO, NLO, NNLO and b) (lower) LL, NLL, NNLL accuracy. 

As expected from the QCD running of as, the cross sections typically decrease when 
increases around the characteristic hard scale Mh, at fixed \ip — Mh- In the case of variations 
of ix p at fixed /xr = M H , we observe the opposite behaviour. In fact, when M H = 115 GeV, 
the cross sections are mainly sensitive to partons with momentum fraction x ~ 10~ 2 , and in 
this x-range scaling violation of the parton densities is (moderately) positive. Varying the two 
scales simultaneously (/xf = /Ar) leads to a compensation of the two different behaviours. As a 
result, the scale dependence is mostly driven by the renormalization scale, because the lowest-order 
contribution to the process is proportional to ccg, a (relatively) high power of a$. 

Figure IBk shows that the scale dependence is reduced when higher-order corrections are in- 
cluded. When resummation effects are implemented (Fig. Eb), we typically observe a further 
(slight) reduction of the scale dependence, with the exception of the factorization-scale depen- 
dence at fixed fi r = Mh that is marginally stronger after resummation. This suggests that the 
rather flat dependence on fxp at NNLO can be an accidental effect, as also suggested by the fact 
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Figure 7: Scale dependence of the Higgs production cross section at the LHC for Mh = 400 GeV 
at a) (upper) LO, NLO, NNLO and b) (lower) LL, NLL, NNLL accuracy. 

that the \ip dependence is much weaker than the dependence at each fixed order (LO, NLO, 
NNLO). 

In Fig. analogous results are plotted at a higher value, M H = 400 GeV, of the Higgs boson 
mass. The overall features of Figs. El and [7| are similar, although we notice that the improvement 
in the scale dependence when higher-order contributions are included is slightly better in Fig. [7| 
than in Fig. |U1 An interesting difference between these two figures regards the [ip dependence 
at fixed fiR. The LO, NLO, NNLO and LL results in Fig. show that the corresponding cross 
sections (very) slightly decrease as \ip increases around Mh- This is because, when increasing 
Mh from Mh = 115 GeV to Mh = 400 GeV, the cross section is sensitive to partons with higher 
values of the momentum fraction x, so that scaling violation of the parton densities can become 
slightly negative. The fact that the parton densities are evaluated in an x-range where scaling 
violation changes sign is also suggested by the change in the slope of the \ip dependence when 
going from LL to NLL and NNLL order. 
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Figure 8: Fixed-order and resummed K-factors for Higgs production at the LHC. 



The impact of higher-order corrections is sometimes presented through the -fC-factors, defined 
as the ratio of the cross section evaluated at each corresponding order over the LO result. The 
.ff -factors are shown in Fig- IH1 where the bands are obtained, as in Sect. 14.11 by varying the scales 
jiR and (simultaneously and independently) in the range 0.5M# < /xf, A*h — 2M#, with the 
constraint 0.5 < Hf/^r < 2. The LO result that normalizes the ^-factors is computed at the 
default scale Mh in all cases. We see that the effect of the higher-order corrections increases with 
Mh- We also see that the soft-gluon resummation effects are more important at higher values of 
Mh. This is expected, since by increasing Mh we are closer to the hadronic threshold, where soft- 
gluon effects are larger. When Mh increases, the scale dependence after resummation is smaller 
than at the corresponding fixed orders. In the case of a light Higgs boson (M#^200 GeV), the 
NNLO K-factor is about 2.1-2.2, which corresponds to an increase of about 20% with respect to 
the NLO K-factor. In this low-mass range, the effects of resummation are also moderate: at NNLL 
accuracy the central value of the cross section increases by about 6% with respect to NNLO. 

In Fig.Elwe plot the NNLO and NNLL cross sections, with the corresponding scale-dependence 
bands (computed as in Fig. |HJ), in the range Mh =100-300 GeV. The corresponding numerical 
results are given in Table [TJ where <T min , er max and <r re f correspond to the minimum, maximum and 
central values in the bands. 



5.3 Tevatron 

Here we study the phenomenological impact of soft-gluon resummation on the production of the 
SM Higgs boson at the Tevatron Run II. 

As in the previous subsection, we show in Fig. El the scale dependence of the fixed-order and 
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Table 1: NNLO and NNLL cross sections (in pb) at the LHC, using MRST2002 parton 
densities. 
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Figure 9: NNLL and NNLO cross sections at the LHC, using MRST2002 parton densities. 

resummed results. We use M H = 150 GeV. As in the LHC case, the cross sections typically 
decrease when //r increases around the characteristic hard scale Mh ■ Figure ITUk shows that the 
fixed-order cross section decreases when fip increases at fixed //r. This is not unexpected: at the 
Tevatron, the cross section is mainly sensitive to partons with x ~ 0.05-0.1, where the scaling 
violation is slightly negative. As in the LHC case, the [ip dependence of the resummed results 
appears to be stronger than the [ip dependence of the fixed-order results. The slope of the /ip 
dependence changes sign in going from LL to NLL and NNLL order, as in the case of Fig. [7| 

In Fig. ^2 we plot the K- factor bands defined as in Fig. |HJ We see that the impact of higher- 
order corrections at fixed Mh is larger at the Tevatron than at the LHC. This is not unexpected: 
at the Tevatron, the Higgs boson is produced closer to the hadronic threshold and soft-gluon effects 
are therefore more sizeable. The NNLO K-factor is about 3-3.2, with an increase of about 40% 
with respect to NLO. Correspondingly, also the resummation effects are larger: at NNLL they 
increase the NNLO cross section by about 12-15%. The scale dependence of the NNLL result is 
smaller than at NNLO. 

In Fig. ^1 we plot the NNLO and NNLL cross-section bands, obtained as in Figs. IH1 and II 11 in 
the mass range Mh = 100-200 GeV. The corresponding numerical results are given in Table El 

5.4 Uncertainties on Higgs production cross section 

In this subsection we would like to discuss the various sources of QCD uncertainty that still 
affect the Higgs production cross section, focusing on the low-M# region (M#^200 GeV). The 
uncertainty basically has two origins: the one coming from still unknown perturbative QCD 
contributions to the coefficient function in Eq. (^Q), and the one originating from our limited 
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Table 2: NNLO and NNLL cross sections (in pb) at the Tevatron (yfs = 1.96 TeV), 
using MRST2002 parton densities. 
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Figure 10: Scale dependence of the Higgs production cross section at the Tevatron for Mh = 150 
GeV at a) (upper) LO, NLO, NNLO and b) (lower) LL, NLL, NNLL accuracy. 

knowledge of the parton distributions. 

Uncalculated higher-order QCD radiative corrections are the most important source of un- 
certainty on the coefficients G a b- A method, which is customarily used in perturbative QCD 
calculations, to estimate their size is to vary the renormalization and factorization scales around 
the hard scale Mh- In general, this procedure can only give a lower limit on the 'true' uncertainty. 
This is well demonstrated by Figs. |H1 and ^2 which show no overlap between the LO and NLO 
(or, LL and NLL) bands. However, the NLO and NNLO bands and, also, the NNLO and NNLL 
bands do overlap. Furthermore, the central value of the NNLL bands lies inside the corresponding 
NNLO bands. This gives us confidence in using scale variations to estimate the uncertainty at 
NNLO and at NNLL order. 

Performing scale variations as in Figs. |H1 El ^2 an d from the numerical results in Tables ^ 
and 121 we find the following results. At the LHC, the NNLO scale dependence ranges from about 
±10% when M H = 120 GeV, to about ±9% when M H = 200 GeV. At NNLL order, it is about 
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Figure 11: Fixed order and resummed K-f actors for Higgs production at the Tevatron Run II 
(V« = 1-96 TeV). 

±8% when M H & 200 GeV. At the Tevatron, when M H & 200 GeV, the NNLO scale dependence 
is about ±13%, whereas the NNLL scale dependence is about ±8%. 

Another method to estimate the size of higher-order corrections is to compare the results at 
the highest available order with those at the previous order. Considering the differences between 
the NNLO and NNLL cross sections, we obtain results that are consistent with the uncertainty 
estimated from scale variations. 

We have also considered other possible sources of higher-order uncertainty. The NNLL coeffi- 
cient is not yet exactly known (see Eq. (|38|) and the accompanying comment). One can thus 
wonder which is the numerical impact of A^ 3 ' in the calculation. We have investigated the nu- 
merical effect of A^ by comparing the full NNLL result with the same result obtained by setting 
A^ = 0. The differences are below 1%, allowing us to conclude that the uncertainty from A^ 
can safely be neglected. Similar conclusions can be drawn about the inclusion of the dominant 
collinear contributions (see Eq. ()48j) ) in the resummed formula. At NNLL order, it gives an effect 
below 1%, in agreement with the fact that collinear logarithmic contributions, which are formally 
suppressed by powers of 1/N, give only a small correction to the dominant terms due to soft and 
virtual contributions (see Sect.EJ). 

A different and relevant source of perturbative QCD uncertainty comes from the use of the 
large- M t approximation in the computation of G a b beyond the LO. The comparison |H1EI] between 
the exact NLO cross section and the one obtained in the large- M t approximation (but rescaled with 
the full Born result, including its exact dependence on M t and M&) shows that the approximation 
works well also for M#^M 4 . This is not accidental. In fact, the higher-order contributions to 
the cross section are dominated by soft radiation, which is weakly sensitive to the mass of the 
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Figure 12: NNLL and NNLO cross sections at the Tevatron (y/s = 1.96 TeV), using MRST2002 
parton densities. 



heavy quark in the loop at the Born level. In other words, as for the size of the QCD radiative 
corrections, what matters is that the heavy-quark mass, M t is actually larger than the soft-gluon 
scale, Mh/N, rather than the Higgs boson scale Mh- This feature, i.e. the dominance of soft-gluon 
effects, persists at NNLO (see Sect. EJ) and it is thus natural to assume that, having normalized our 
cross sections with the exact Born result, the uncertainty ensuing from the large- M t approximation 
should be of order of few per cent for Mh ^200 GeV, as it is at NLO. 

Besides QCD radiative corrections, electroweak corrections also have to be considered. For a 
light Higgs, the 0{GfM^) dominant corrections in the large-Mt limit have been computed and 
found to give a very small effect (well below 1%) [52] . 

The other independent and important source of theoretical uncertainty in the cross section is 
the one coming from parton distributions. 

We start our discussion by considering NLO results. We have compared the MRST2002 NLO 
cross sections with the ones computed with CTEQ6 [53] and Alekhin [53] distributions. All three 
sets include a study of the effect of the experimental uncertainties in the extraction of the parton 
densities. At the LHC, we find that the CTEQ6M results are slightly larger than the MRST2002 
ones, the differences decreasing from about 2% at Mh = 100 GeV to below 1% at Mh = 200 GeV. 
Alekhin's results are instead slightly smaller than the MSRT2002 ones, the difference being below 
3% for Mh ^ 200 GeV. At the Tevatron, CTEQ6 (Alekhin) cross sections are smaller than the 
MRST2002 ones, the differences increasing from 6% (7%) to 10% (9%) when Mjj increases from 
100 GeV to 200 GeV. These results are not unexpected, since in Higgs production at the Tevatron 
the gluon distribution is probed at larger values of x, where its (experimental) uncertainty is 
definitely larger. The cross section differences that we find at NLO are compatible with the 
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Table 3: NNLO and NNLL cross sections (in pb) at the LHC (fip = [Ir = M^) using 
the parton distributions of Ref. \54Jj - 



experimental uncertainty on the NLO gluon luminosity quoted by the three groups, which is 
below about ±5% at the LHC and about ±10% at the Tevatron |_3ES1IHI]- 

Throughout the paper we have used the MRST2002 set as the reference set of parton 
distributions. The main motivation for this choice is that this set includes (approximated) NNLO 
densities, allowing a consistent study at NNLL (NNLO) accuracy. The CTEQ collaboration does 
not provide a NNLO set, and a complete, consistent comparison with MRST can therefore not 
be performed. In Ref. a NNLO set of parton distributions (set A02 from here on) has been 
released. This set also includes an estimate of the corresponding uncertainties. We can thus 
perform a comparison with the MRST2002 results. 

In Table _] (Table _J we report the NNLO and NNLL cross sections (computed with fiF = 
fifi ; = Mh) obtained at the LHC (Tevatron) by using the set A02 of parton distributions. We 
also report the corresponding errors, computed by using the dispersion array* of the A02 sets. 
Comparing Tables _] and _] with the corresponding central results in Tables _] and _J we see that 
there are relatively large differences that cannot be accounted for by the errors provided in the 
A02 set. At the LHC the cross section is larger than the one computed with the MRST2002 set, 
and the difference goes from about 8% at low masses to about 2% at Mh = 200 GeV. At the 
Tevatron the cross section is lower than that using MRST2002, with a difference that ranges from 
about 7% at low M H to about 14% at M H = 200 GeV. 

These differences can be better understood by looking at Figs. __] and __J where a comparison 
of Alekhin (with errors) and MSRT2002 luminosities is presented. Comparing Tables _J _J _] and _] 

*More precisely, 30 parton distributions are generated according to fi(x, k) = fi(x) ± A^(x, k), k = 1, . . . , 15, 
and the cross section for each k is evaluated. The positive (negative) deviations from the central value are then 
summed in quadrature to obtain the upper (lower) error. 
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Table 4: NNLO and NNLL cross sections (in pb) at the Tevatron ([ip = fiR = Mjj) 
using the parton distributions of Ref. J5^| /. 



with Figs.EHandEl we see that the differences in the cross sections are basically due to differences 
in the NNLO gluon-gluon luminosity. Moreover, Figs. ITBI and IT4l show that the differences between 
the luminosities are typically larger than the estimated uncertainty of experimental origin. In 
particular, the differences between the gg luminosities appear to increase with the perturbative 
order (i.e. going from LO to NLO and to NNLO). 

We are not able to pin down the origin of these differences. References |5U] and [SU use the 
same (though approximated) NNLO evolution kernels, but the MRST2002 set is based on a fit 
of deep- inelastic scattering (DIS), DY and Tevatron jet data, whereas the A02 set is obtained 
through a fit to DIS data only. 

From the above discussion, we conclude that the theoretical uncertainties of perturbative origin 
in the calculation of the Higgs production cross section, after inclusion of both NNLO corrections 
and soft-gluon resummation at the NNLL level, are below 10% in the low-mass range. However, 
it is also apparent that there are uncertainties in the (available) parton densities alone that can 
reach values larger than 10%, and that are not fully understood at the moment. Improvements 
of these aspects may only come from better understanding of parton density determinations, and 
are therefore totally unrelated to the calculation presented in this work. 



6 Conclusions 



In this work we have presented a next-to-next-to-leading resummation of soft-gluon effects in 
the gluon fusion cross section for Higgs production. Furthermore, we have supplemented our 
resummed result with available fixed-order (up to NNLO) calculations. We have presented a 
phenomenological study of the Higgs cross sections at the LHC and the Tevatron, and showed 
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Figure 13: Comparison of A02 and MRST2002 gg luminosities at the LHC. 

that, once problems with parton density determinations are overcome, a theoretical precision of 
about 10% in the predicted cross section can be attained. 

The impact of the NNLL corrections included in this work was found to be modest in the case 
of light-Higgs production at the LHC. The NNLL corrections are more significant at the Tevatron, 
and allow us to reduce the theoretical error band at the same precision level as at the LHC. We 
believe that our results add much confidence in the predicted cross sections for the following 
reason. The pattern of radiative corrections from the NNLL resummation, when truncated at 
fixed order, reproduces quite well the pattern from the exact fixed-order calculations. This adds 
confidence that the higher-order (beyond the NNLO level) radiative corrections arising from the 
NNLL treatment reflect the behaviour of the exact theory. 

As a side product of this work, we have shown that the iV-space formulation of the soft-gluon 
approximation is superior to the x-space one. In Ref. j^H] it was shown that the iV-space formu- 
lation of soft-gluon resummation should be used, since it avoids a gross violation of momentum 
conservation, which in turn leads to an unphysically large asymptotic growth of the coefficients of 
the perturbative expansion. In the present work, the better behaviour of the iV-space approach 
has also been shown to hold in practice, when comparing known fixed-order results with the 
fixed-order expansion of the resummation formulae, up to the NNLO level. 



A Appendix: iV-moments of soft-gluon contributions 

In this appendix we present a simple method to evaluate iV-moments of soft-gluon contributions 
to arbitrary logarithmic accuracy. 
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Figure 14: Comparison of A02 and MRST2002 gg luminosities at the Tevatron. 



We consider the iV-moments I n (N) of the singular distributions T> n (z) defined in Eq. (fTTj) : 



r 1 z N ~ l - 1 



I n {N) = I dz z N ~ 1 V n (z) = I dz ln n (l - z) . (63) 



These iV-moments were computed in Ref. jzo]. It was shown that, in the large- N limit, the leading 
and next-to-leading logarithmic terms (ln n+1 N and In™ N) arising from the integration over z can 
straightforwardly be obtained by a very simple prescription. It is indeed sufficient to replace 
the integrand weight (z 1 ^" 1 — 1) by the approximation 

z*- 1 - 1 ~ - 6(1 - * - JVo/iV) , (64) 

where A^o = e~ lE . In the following we show how this NLL prescription can be generalized to any 
logarithmic accuracy. The all-order generalization of Eq. (|E4*|) is given in Eq. (f?4*|) or, equivalently, 
in Eq. (fTHJl. 

The iV-moments in Eq. (JBlIj) can be evaluated as described in Ref. Using the following 
identity 

d 



\n n (l-z) = hm( — ) (1-zY (65) 



to replace the logarithmic term in the integrand on the right-hand side of Eq. (|63J). we obtain 

T(N) r(l + e) 



(66) 



where T(x) is the Euler T-function. The expression (fMJ) can be used to evaluate I n {N) for any 
given value of n without approximations of the dependence on iV 
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We are interested in the large- N behaviour of I n (N). Neglecting terms that are suppressed by 
powers of 1/N in the large- N limit, we can use the approximation 



r(N + e 

and we can rewrite Eq. (Jfifij) as 



r W _ e -,ln« 



(67) 



= S- (|)" {i[«—r(i + e) -i] 



Using the expression 



r(l + e)=exp|- TE e + £(-ir^e"| , (69) 



the term in the curly bracket of Eq. (J68j) can easily be expanded in powers of e. The result for 
I n (N) is thus a polynomial of degree n + 1 in the large logarithm lniV: 

I n (N) = [ -^— (InTV + lE ) n+l + n C(2) (In JV + t^)™" 1 

71 + 1 2 

n-2 

(70) 



n-2 , , 

+ J2 *nk (lnN + lE ) k + 0(-) 

k=0 ^ ' 



where the coefficients a n k are combinations of powers of C(2), C(3), • • • , C( n + 1) ( see e -g- Table 1 
in Ref . 



As shown by Eq. (J70D and pointed out in Ref. [H], the LL and NLL terms (ln n+1 N and ln n N) 
can straightforwardly be obtained from the definition of the iV-moments in Eq. (|63|) by simply 
implementing the prescription in Eq. (|64[). To show how this NLL prescription can be generalized 
to any logarithmic accuracy, we consider the right-hand side of Eq. ()68|) and use the following 
formal identity: 

e~ €lnN T(l + e) = rf 1 - — — } e~ elnN . (71) 



dm TV, 

Then we can perform the n-th derivative with respect to e, and we obtain 



This expression can be regarded as a replacement of Eq. (joHjl to compute the polynomial coefficients 
a n k in Eq. (|70|l . Moreover, by noting that 



71+1 Jo 1 — 2 

we shortly get the all-order generalization of Eq. (|64j). It is: 

z^- 1 - 1 = - r(l - 9(1 - z - 1/JV) + C(1/JV) (74) 

= - f(l - ^y) 6(1 - z - iVo/iV) + 0(1/JV) , (75) 
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where we have also introduced the function T defined by 



+00 



r(l + e) = e^ e r(l + e)=exp|^(-l) n ^e"| . (76) 

Note that the equality between Eqs. (|71jl and (|75j) simply follows from T = T exp [y E Nd/(dN)} 
and from the fact that the exponential of the differential operator acts as the rescaling iV — > 
exp[ lE }N = N/Nq. 

Although we have derived it by starting from Eq. (|63|). it is straightforward to show that the 
prescription ((74")) (or ((75))) can be applied as follows 

1 z n-i _ 1 / a \ /-l-i/JV . / j 

rfz __ nas , ln(1 _, ))= _r(i-_) I _ F(as , ln(1 _ z)) + (_ 

to evaluate the In TV-contributions arising from the integration of any soft-gluon function F that 
has a generic perturbative expansion of the type 

+00 2k— 1 

F(t* s , Ml - ^)) = E E ^"(1 - z ) • ( 7? ) 

fc=l n=0 



The all-order differential operator T (T) in Eq. (|74|) (Eq. ()75|)) is obviously defined only in a 
formal sense through Eq. (J69|) (Eq. (J76|) ) . However, such a definition is sufficient for any formal 
manipulation to all logarithmic orders. Moreover, it is also sufficient for any practical use at fixed 
and arbitrary logarithmic accuracy. Since each additional power of d/(d\nN) leads to terms that 
are suppressed by a power of IniV, we can obtain all the terms up to a given (and arbitrary) 
logarithmic order by simply truncating the power series expansion of T at the corresponding order 
in the logarithmic derivative. For example, using the truncation 

+ i(7i + 3 7 EC(2) + 2C(3))(^)V.. , (78) 

the first, second, third and fourth terms on the right-hand side lead to the LL, NLL, NNLL and 
N 3 LL contributions, respectively. 



B Appendix: Equivalence between resummation formulae 



Here we prove the equivalence between the two all-order representations (J29|) and (|34j) of the 
Sudakov radiative factor A^. In particular, we derive the relations between the functions A(as) 
and -D(as) in Eq. (j2S]) and the functions D(as) and C gg (as) in Eq. J33J). These relations are: 



D(a s ) = D(a s )+4d as T 2 (d, 



as , 



A(a s ) - - d 



D(a s ) 



(79) 
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C M (a s (M|),M|/^ = 1) = exp |-4 T 2 (d t 
where the differential operator d as is defined by 



A(a s ) - j d as D(a s ) 



(80) 



d as = -2p(a s ) a s 



/3(as) is the QCD /3-function, 



d 



— — — = P{a s {n )) = -}_^b n a s , 



32) 



n=0 



whose first three coefficients, bo,bi,b2, are given in 
its power series expansion in e: 



r*(e) = I [i 



-7 B e 



r(l-e)] = -H-exp 



Eq. (jlHJ), and the function T 2 (e) is defined by 

(83) 



boo 



1=2 



Note that Eq. JSDJ) gives C gg (as(n 2 R ), M 2 H j y? R ) at the scale /i| = M\ as a function of a s (M^). The 
dependence on fj, R is simply recovered by using renormalization-group invariance, i.e. by expressing 
as{Mff) as a function of as{fi 2 R ) and Mfj/fi 2 R through the solution of the renormalization-group 
equation (jH2) (see also Eq. Pty). 



Inserting the following expansion 

r 2 (6) = -k(2)-^C(3)6 + 0(e 2 ) , 



in Eqs. (179)1 and (jHUj) . we can straigi 
up to NLL accuracy: 



2 ^ v ' 3 

;htforwardly obtain the expansions of the functions D and C 



14) 



99 



5(a s ) - 2?(a s ) = - (^) 2 4C(2)ttM (1) + 0(a 3 s ] 
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^(«sK),M^)=l + 



if' 



7T 



2C(2)A« 



a s (/if ^ 



OR j 



71 



2C(2) ^ 2 )-7rM (1) ln^) 



+2 (C(2)A«) 2 + ^C(3) 7 rM {1 



(86) 



To derive Eqs. (|75j) and ()80|) . we start from Eq. ()29)1. we replace the integrand weight (z N ~ x — 1) 
by the all-order logarithmic prescription in Eq. 1)75)1. and we change the integration variable z as 
^ — > y — 1 — z. We obtain: 



In A* os(a*^ 



' 2 ' 2 



r i 



d 



dy 

dlnNj J No/N y 



y 2 M 2 H 



I-'f 



dq 2 

.2 ]i j2 



, >%s(<7 2 )) (87) 
D{a s {y 2 M' 2 H )) +0(1/N) . 
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Then we use the definitions in Eqs. (J75|) and (jHHJ) to get 







\ dhiN 



i-r, 



o 



d\nN \d\nN 







Inserting Eq. ()88)1 in Eq. (jHTjl and comparing the latter with Eq. (|34jl . we get 



Vo/AT 2/ 

-r 2 



D(a s (y Mfj)) — \nC gg (a s (Mf I ), Mfj/^ = 1 



nq/n y 



D(a s (y 2 M 2 H )) 







dhiN \d\iiN 







2 r x 



dy 

N /N y 



dq 2 



A(a s (q 2 )) +D(a s (y 2 Mf I )) 



19) 



f ^ D(a s (y 2 M 2 H )) 

jno/n y 

-4A(a s (N 2 M 2 H /N 2 )) 







dlniV 







dhiN 



D(a s (N 2 M 2 /N 2 )) 



(90) 



where Eq. ()90)) is obtained from Eq. (|5U|) by acting with the operator ( g^^ ) 2 onto the y-integral. 
Using the renormalization-group equation (|82l). we now observe that the relation 



5 



<91niV 



/(os(fcV^)) 



5 



(9 In & s 



/(as(fcV^)) = 9 Qs /(a s (fcV^)) 



(91) 



is valid for any arbitrary function /(as) (the operator <9 as is defined in Eq. (|81|)). Therefore, we 



can perform the replacement 



a In AT 



d as in Eq. (j9*U)l . and we obtain 



dy 

N a /N V 



D(a s (y 2 M 2 H )) - In C gg (a s (M 2 H ), M 2 H /pi 2 R = 1 



ATq/AT 2/ 



D(a s (y 2 M^)) 



r 2 (a as ) 4A(a s ) - <9 QSJ D(« S ) 



«s=as(A r 2M lr/ Ar2 ) 



(92) 



Note that this equation holds for any value of N. Thus, setting N = No the y-integrals vanish, 
and we immediately get Eq. (|8U|) . Moreover, we can apply the operator dl ® N on both sides of 

Eq. fj92jl and, since C gg is N- independent, we obtain the relation 



D(a s (N 2 M 2 H /N 2 )) = D(a s (N 2 M 2 H /N 2 )) 

d 



+ 



d\nN 

which, owing to Eq. (|9~TJ). gives exactly Eq. (|79*|) 



r 2 (a QS ) 4A(« s )-a asJ D(a s ) 



a s =as(NgMf I /N 2 ) 



(93) 
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C Appendix: Resummation formulae at NNLL accuracy 
and beyond 



In this appendix we sketch the calculation of the logarithmic functions g^ in Eq. ([27)1 . These 
functions originate from the Sudakov radiative factor in Eq. ()28|) . Thus we write 



lnA^f as (A* 



4), 



M H. M H 



f4 



InN g£\b a s ^ R )lnN) + 



(94) 



+ E [«s(4)] n " 2 gf{b a a^l)\nN,M%/^Ml/^ F ) + 0(1) 



n=2 



where the term 0(1) stands for contributions that are constant in the large-iV limit. This TV- 
independent term eventually contributes to the function C gg (as) in Eq. 



The logarithmic expansion on the right-hand side of Eq. (|9~4"j) is most easily computed by 
starting from the integral representation in Eq. (J34j) and by using the method described in Refs. [2H1 
IHH] . We first use the renormalization group equation (JH2J) to change the integration variables {q 2 , y} 
in {as(q 2 ) = a,as(y 2 M]j) = a 1 }, so Eq. (pMjl becomes 



InAf 



Ml 

/4 



<xs(M%) 



da' 1 



a s (N$Ml/N2) «' Pi®') 



da A{a) 



<*s(/4) a 



D(a) 



O(l). 

(95) 



Now, the integrand can be expanded (see Eqs. (jSHJ), (JHH), (JZHJ) and (JB2J)) in power series of a and 
a', and the expansion can be truncated to the required (and arbitrary) logarithmic accuracy. This 
procedure leads to simple (logarithmic and polynomial) integrals, and the result is expressed in 
terms of elementary functions of the perturbative coefficients A^ n \ D^ n \ b n and of a$(k 2 ) with k 2 = 
NqM]j/N 2 , M 2 h or fj,%. To obtain the functions g^' '(^asO^) IniV, Af£//4; Mj^//4), it is sufficient 
to express as{k 2 ) in terms of as{fi 2 R ) and ln(/c 2 //^) (i.e. lniV,ln(M#//x 2 ? ),ln(Mjj//4i),ln iVo = 
—Je) according to the perturbative solution of the renormalization group equation (J82|) . and then 
to compare the results with the right-hand side of Eq. ()94j) . For instance, to extract the LL, NLL 
and NNLL functions g$,gP and g^§ in Eqs. pi jl .p ty and (jUJ, it is sufficient to use the NNLO 
solution of Eq. 



«s(/4) 



bo 



0(al^ 2 R )(as(^ R )Hk 2 /^ R )r 



b 2 
°o 



-ln£ + 



\ bo 

1 _ r 



(96) 



where we have defined £ = 1 
accuracy is straightforward. 



■bo as(fj, 2 R ) ln(k 2 / fi 2 R ). The extension to arbitrarily higher logarithmic 
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n, 2 

a s 


a s 


a s 


a s 


a s 


n> >7 


LHC 

Tevatron 


11.67 
0.228 


14.72 
0.291 


8.61 
0.194 


1.65 
0.065 


0.30 
0.017 


0.08 
0.007 



Table 5: Contributions to the total cross sections (in pb) at the LHC and the Tevatron 
from higher orders in the expansion of the NNLL resummed result, with [ir = fip = Mh , 
M H = 130 GeV and MRST2002 NNLO parton densities. The sum of the first three 
columns gives the exact NNLO result. 



D Appendix: Convergence of the fixed-order 
soft-gluon expansion 

We check here the numerical convergence of the fixed-order soft-gluon expansion toward the re- 
summed result obtained by using the Minimal Prescription. 

The rapid convergence of the higher-order soft-gluon corrections is displayed in Table for 
the case of Higgs production at the LHC and the Tevatron. All the entries reported in Table 03 
are obtained by using the MRST2002 NNLO parton densities. The first three columns show the 
contributions to the total cross section from the first three orders in perturbation theory. The 
sum of the first three columns thus gives the exact NNLO result. The next two columns show 
the contributions from the next two fixed-order terms, as obtained by the (truncated) fixed-order 
expansion of our NNLL resummed formula. The sum of all entries in each row corresponds to the 
full NNLL resummed result as obtained with the Minimal Prescription, while each fixed-order term 
has no ambiguity due to the choice of the contour for the Mellin transformation in Eq. (|61|). The 
results in Table El support the validity of the Minimal Prescription, since the truncated resummed 
expansion converges to it very rapidly. 
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E Appendix: Soft-gluon expansion at N 3 LO 



In this final appendix we present the soft-gluon approximation, ^ N , in N space of the N 3 LO 



(3) 



, ^gg^Ni t° the coefficient function in Eq. (0). The analytic expression of G N 



contribution, G 

obtained by performing the expansion of the all-order resummation formula in Eq. (|24|) . is 

G (3)S V-iV (M 2^| )M 2 //4) = H^pl ln * N 

- \8 (A^) 2 b 7T 



(3) SV-AT 



+ 8lE (A«) 3 -4 (A«) 3 ln^ 

Lip 



+ ln 4 A7 



4^(1)^(2) 



4A^b 2 TT 2 



+ 2(AUyC% + ±im ^ + 16(^))M 



- 4 (^)) 2 ^ In # + 4 (^)) 3 In 2 ^ - In ^ ( 8 J^1 + 16 (A (D) 
A*b A*f A*f V ^ 



+ In A/" 



3 



A*F 

(1) 



7e 



8^ + + «A<%,CJ _ 2 c(2) + g {Am y^ c(2) 



+ 



32 (A«) 3 7 i 



24 (^yfeoTTTl + TF^ie^ 1 )^ 2 ) 



16 AW 6 2 



7T 



4 



In 3 —g- + In 2 



A*F 



■ 3 2 

7f 



Mf 
Mj 

2 



ln M|^8AWfe 2 7r 2 
/4 



+ 16 (^ 1 )) 2 6 vr 7s )-ln^(8A«A( 2 )+4 (A«) 2 C« + 8 (A«) 2 b ir lE 
+ 16 



/ , , Ml , Ml 
8 6 tt ln-f ln-f 

Lip Ll R 



+ \n 2 N 



2 A< 3 > + 2 A( 2 > C« + 2 Cg) - 2 b 7r + 8 6 2 vr 2 C(2) + 16 (A (1) ) 2 6 vr 7 | 
+ 7 I (8 (A^) 2 C« + 16 ^ 2 > + 8 6 2 vr 2 ) + 7 f f 8 A (2) 6 tt + 4 A« M 2 



+ 4A«6 7rC« - 4 A<»D<® + 16 (A (1) ) 2 ^C(2)) + 2 A« 6 2 tt 2 In 



.2^1 
/4 



zMl 



-4 (AW) 2 M ln 2 ^ ln^ + 16(A«) 2 b 7r lE ln^f + 2 {A^f b Q n In 

A*F A*B A*F A*B ^F 

+ \n 2 ^f(4A^A^ + 2(A^) 2 C^-4 (A^) 2 b 7r lE ) + In M ( 2 A« D< 2 > 



/^F 



-16A«A( 2 ) 7i5 -8 C« 7 f-8 (A«)^ 7r 7 I-8 

- In ^ (4 A< 3 > 6 tt + 2 A« 6 l7 r 2 + 2 b n C« + 8 6 2 tt 2 7 f + 16 (A« ) 2 b 7r] 
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+ lniV 



- - Cg) + 32A{1) b l!llM + 8 6o7r c(2 ) + 4 M 2 C(2) 

3 

+ 4 ^ &ott C$ C(2) + 1QA(1) ^ 2 1 3 e + 7 2 ^ A ( 2) ^ + 4 A (x) ^ + 4 &( 
+ 7S (4 A® + 4 A< 2 > C« + 4 AW Cg) - 4 6 vr £>< 2 ) + 16 A™ 6 2 vr 2 C(2) 



2A« 6 2 



lir 



M 



+ 2 A« & 2 tt 2 In 2 ^ In ¥f - 2 6 2 vr 2 In In 2 ^ + 4 v4 (1) 6 2 vr 2 7£ In 2 ^ 

A*F A*i? A*F f^R f^R 

_ i n 2 ^ ( 2 A m bQ7T + A (D bin 2 + A (i) ^ ^(lA _ i n ^ ( 2 A< 3 > + 2 A( 2 ) C« + 2 Cg)) 



+ ln Mf ln I lfU A® bo7T + 2 A« M 2 + 2 A (1 ) M C« 
Mi Mi V 99 



M 



/4 



/4 



- In ^ (s A( 2 ) & tt 7s + 4 v4 (1) 6 l7 r 2 7s + 4 v4 (1) C« & tt 7 £ + 8A^ 6 2 vr 2 T | + S^ 1 ) 6 2 tt 2 C(2) 

f^R V 



(97) 



Note that with the present knowledge (see Sect. l3.2|) of the perturbative coefficients of the functions 
A(as), D(as) and C gg (as), it is possible to predict the numerical values of the coefficients in the 
expansion (J97j) only up to 0(ln 2 N) terms. The term proportional to In N and the ^-independent 
term depend on the unknown coefficients and Cf g \ respectively^. The same results apply 
to the DY process, by simply replacing the gg-channel coefficients A^ n \D^ and with the 
corresponding gg-channel coefficients. The expression (}9~Tj) can be used to check future perturbative 
calculations at N 3 LO. 
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